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1.  Introduction 


Brittle  materials  (i.e.,  ceramics)  are  being  considered  for  an  increasing  number  of  applications 
ranging  from  armor  elements  to  gun  tube  inserts  to  engine  components,  where  they  must 
withstand  intense  dynamic  loads.  The  damage  and  massive  failure  processes  that  can  occur  in 
these  conditions  are  not  well  understood.  Because  ceramics  typically  have  lower  symmetry  than 
most  metals,  the  basic  crystal  components  are  not  only  highly  anisotropic  but  also  have  only 
limited  slip  systems  for  plastic  deformation  in  polycrystalline  configurations.  Therefore,  they 
must  rely  on  other  mechanisms  such  as  twinning,  microcracking,  and  void  formation,  as  well  as 
limited  plasticity,  to  achieve  required  deformation.  Adjacent  grains  of  differing  orientation  will 
find  it  hard  to  maintain  coherent  boundaries;  large  stress  concentrations  within  grains  may  build 
up;  and  large  oscillations  in  stress  and  strain  may  occur  from  grain  to  grain.  Homogenization 
theories  can  predict  average  material  properties.  However,  in  the  circumstances  just  described, 
knowledge  of  the  oscillations  is  also  necessary.  It  is  the  purpose  of  this  effort  to  develop 
techniques  that  will  give  an  account  of  the  oscillations,  as  well  as  the  averages,  and  to  apply  the 
resulting  theory  to  several  important  topics  such  as  damage,  incipient  microcracking,  and 
fracture  waves. 

During  several  of  the  last  decades,  ceramic-based  annor  technologies  and  corresponding 
experiments  and  theories  have  shown  a  remarkable  progress.  One  of  the  brightest  achievements 
of  the  last  decade  is  the  discovery  of  the  phenomena  of  dwell-penetration  transition  and  of  a 
defeat  of  a  projectile  targeting  ceramic  armor  (1-3).  The  phenomenon  of  a  projectile  defeat  is 
accompanied  by  the  appearance  of  an  intensively  fractured  zone  (IFZ)  within  the  bulk  of  a 
ceramic  target.  The  shape  of  this  zone  resembles  a  tri-axial  ellipsoid  (see  figure  1,  where  the  IFZ 
zone  is  shown  in  white).  In  their  recent  publication,  LaSalvia  and  Normandia  (4)  interpreted  it 
on  the  basis  of  the  famous  Bussinesq  solution  of  classical  theory  of  elasticity.  Similar 
observations  of  IFZ  are  known  in  ceramic  indentation  experiments  (5)  and  in  the  seismology  of 
earthquake  epicenters  (6). 

When  dealing  with  IFZ,  it  is  often  desirable  to  extend  the  information  about  the  overall  modulae 
with  the  information  about  the  shape  and  morphology  of  the  IFZ.  Such  information  can  be 
obtained  by  using  additional,  specific  criteria  or  more  general  principles  of  mechanics  or 
thennodynamics.  In  this  report,  following  the  minimum  energy  approaches  of  Gibbs  (7-8)  and 
Griffith  (9),  we  propose  one  such  thennodynamic  approach,  which  allows  one  to  determine  the 
shape  of  the  IFZ  simultaneously  with  the  distributions  of  stresses  and  strains  within  the  systems 
of  interest.  The  general  methods  and  results  of  the  suggested  approach  are  illustrated  by 
considering  an  instructive  problem  of  nucleation  of  a  damaged  zone  in  an  isotropic  solid  space 
under  sufficiently  high  uniform  stresses  at  infinity. 
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Source:  Picture  courtesy  of  J.  LaSalvia  and  M.  Normandia  (4). 


Figure  1.  Damaged  observation  from  recovered  ceramics. 


2.  The  Substance  Model 


Consider  a  macroscopically  heterogeneous  configuration  possessing  two  macrodomains,  co+  and 
co_  ,  occupied  by  the  consolidated  and  damaged  modifications,  respectively.  Let  y  be  the 
interface  separating  the  two  domains  and  S  the  external  boundary  of  the  system  (see  figure  2). 


Figure  2.  The  geometry  of  the  system. 


The  consolidated  (intact)  modification  is  characterized  by  the  elastic  energy  density  per  unit 
volume 


e+(uiy)  = 


UfJU 


k\l  9 


(1) 


where  u/x)  is  the  displacement  at  the  point  xl(uij  =  du/dx1)  and  cf  is  the  elasticity  tensor  of 

the  consolidated  modification.  The  absence  of  the  free-constant  and  linear  terms  in  the  Taylor 
expansion  of  e(uf  j )  indicates  that  the  configuration  with  n,  =  0  is  stress  free,  and  it  is  chosen  as 
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the  base  for  the  calculation  of  the  elastic  energy  of  the  substance.  For  the  damaged  modification 
of  our  system  we  can  use  an  “effective”  elastic  energy  density  as  follows: 


e_{uiJ)  =  e°  + 


u,  u, ... 

i\j  k\l 


(2) 


The  damaged  modification  is  elastically  weaker  than  the  intact  one,  i.e.,  its  elasticities  in  some 
sense  are  smaller  than  the  elasticities  of  the  intact  modification.  Besides,  there  is  a  nonzero 
constant  e°  in  the  elastic  energy  density  of  the  damage  modification — this  constant  reflects  the 
work  required  to  produce  defects  in  the  intact  modification.  Schematically,  the  two  energies  are 
shown  in  figure  3. 


Figure  3.  The  specific  energy  density  of 
the  intact  ( e+  )  and  damaged 

( e  )  modifications. 


The  parabola  of  the  consolidated  substance  passes  through  the  origin  and  has  the  fastest  growth 
with  e  .  It  is  clear  that  for|e|  <  |e*| ,  the  elastic  energy  density  of  the  consolidated  phase  is  less 

than  in  the  damaged  phase  (for  |e|  >  |e*| ,  the  elastic  energy  density  of  consolidated  phase  is 

greater  than  in  the  damaged  phase).  The  trade-off  between  the  opposite  effects  of  the  damage  on 
“free-stress”  energy  e  and  on  the  elasticity  tensor  cykl  at  various  deformations  is  the  central 

conflict  under  study  with  the  help  of  our  suggested  model. 


The  simplest  model  of  the  previously  mentioned  type  is  the  model  of  an  isotropic  two-state 
elastic  substance.  For  this  model,  the  energy  densities  for  two  modifications  are  given  by  the 
following  formula: 


(Ui\i)  =  K++Hl 


1  2v 


u\,u{.  +U/;,M 


\i  \j 


<i\j) 


(3) 


where  ju±  are  the  shear  Lame  modules  and  v±  are  the  Poisson’s  ratios  (  k  =  0). 
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3.  Continuum  Model  of  Fracture 


Per  classification  of  Grinfeld  and  Wright  (10,  11),  the  model  just  considered  is  a  certain  two- 
state  discrete  fracture  model.  A  continuum  fracture  model  should  include  a  damage  parameter  x 
that  can  assume  a  continuum  range  of  values.  Corresponding  elastic  energy  e  function  is 
becoming  a  function  of  the  elastic  gradients  w/|;  and  of  the  damage  parameter  x  :  e  =  e(ujrx)  . 

For  such  a  model,  the  associated  equilibrium  master  system  includes  the  following  two  bulk 
equations: 

8  de<V)-„ 


de(uiU,x) 


The  associated  system  of  quasistatic  evolution  reads  as  follows: 

d  de(um]it ,  k)  _  () 
dxj  du^ 


dx  _  K  de(ui  rK) 
dt  dx 


where  A  is  a  certain  positive  kinetic  constant  or  function. 

The  associated  dynamic  system  reads  as  follows: 

d2u‘  _  d  de(uniln,K) 
^  dt 2  dx1  du;[1 


ok  _  ^  de(ui  f ,  k) 

dt  dx 


The  important  issue  here  is  the  choice  of  an  appropriate  energy  function  e  =  efio  , ,  x)  and  the 

kinetics  function  K.  The  proper  choice  is  a  hard  problem  of  physics  of  ceramic  materials.  Some 
insight  can  be  acquired  by  using  approaches  accepted  in  phenomenological  damage  theory  (12) 
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and  the  theories  of  proliferation  of  vacancies  in  ceramics  sintering  (13).  Needless  to  say, 
cracking  is  not  a  proliferation  of  vacancies  because  of  various  reasons,  but  this  might  be  a  good 
first  step. 


The  quasistatic  system  presented  in  equations  6  and  7  was  used  to  numerically  explore  the 
behavior  of  a  circular  plate  with  a  thin  elliptic  cavity  under  action  of  external  pressure  (figure  4). 
The  kinetics  function  K  was  assumed  a  certain  positive  constant.  For  the  energy  density,  we 
have  chosen  the  following  Kachanov-Lifshitz  function: 


e(uiU,K) 


f 

&>(*) 

V 


v 

l-2v 


vy.j 


(10) 


with  the  damage  function  </>(k)  (10,  11)  as  follows: 


¥K)  =  \-(\-cmm)y, 

/c 

0 <K<K, 


and 


0  <  c,  ■  <  1 . 


(11) 

(12) 

(13) 


Tm 


Figure  4.  Development  of  damage  in  a  circular  elastic  plate  with  a  cavity. 

Thus,  the  suggested  model  depends  on  six  constants — pi,  v,  £  k  ,  k  ,  and  cmm  .  The  first  two  of 
them  are  just  the  shear  modulus  and  Poisson’s  ratio  of  the  intact  substance;  the  physical  meaning 
of  the  remaining  four  constants  are  explained  elsewhere.  For  a  certain  choice  of  the  constants, 
the  evolution  of  the  relative  damage  parameter  k/k  is  shown  in  figure  4.  The  depth  of  gray 
reflects  the  magnitude  of  the  relative  damage  parameter,  white  corresponds  to  a  fully  damaged 
material,  and  black  to  a  fully  undamaged.  Not  surprisingly,  the  largest  damage  occurs  at  the  tips 
of  the  cavity. 
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4.  Summary  and  Conclusions 


•  The  proposed  theoretical  approach  to  brittle  fracture  allows  the  determination  of  the 
morphology  of  damaged  zones  together  with  the  traditional  fields  like  deformation  and 
stress. 

•  Similar  to  the  Griffith  model,  the  proposed  approach  is  based  on  the  trade-off  between  the 
elastic  energy  and  the  energy  of  chemical  bonds. 

•  In  the  proposed  approach,  the  damaged  and  intact  substances  are  treated  as  different  phases 
of  the  same  substance,  and  “phase  transformation”  damage-recovery  is  taken  into  account. 

•  The  equilibrium  system  includes  coupled  equations  for  the  unknown  displacement  field, 
damage  parameter,  and  morphology. 

•  Because  of  its  deep  nonlinearity,  even  for  small  deformations,  the  equilibrium  system  may 
have  multiple  solutions. 

•  The  equilibrium  system  does  not  allow  static  solution  for  some  boundary  loads;  hence,  the 
system  should  be  extended  by  the  inclusion  of  inertia  and  dissipation  mechanisms. 
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